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ESD-TDR-63-649 


INITIAL  ORBIT  DETERMINATION  FROM  LEAST 
SQUARES  REDUCTION  OF  GEOCENTRIC  POSITION  VECTORS 


ABSTRACT 


A  computer  program  for  initial  orbit  determination  based  on  a 
theory  employing  a  linear  least  squares  fit  to  the  geocentric  position 
vectors  has  been  developed.  Observations  can  be  spread  over  a  number  of 
days  since  this  prodecure  does  not  use  the  times  of  observations;  however, 
they  are  used  in  the  differential  correction  procedure.  This  program 
coded  for  the  Philco  2000,  is  designed  for  operational  use  at  the  Space 
Detection  and  Tracking  Center  in  Colorado  Springs,  Colorado. 
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SECTION  1 


INTRODUCTION 


The  computer  program  discussed  in  this  report  computes  and 
differentially  corrects  the  preliminary  orbital  elements  of  a  satellite 
when  given  an  aggregate  of  observation  position  fixes  (maximum  of  150). 

The  Initial  Orbit  Least  Squares  Program  (IOLSQ)  is  coded  for  the  Philco 
2000  computer  and  is  intended  to  be  a  part  of  the  Semi-Automatic  Program 
System  (SPS)  at  the  SPADATS  Center  in  Colorado  Springs,  Colorado.  This 
program  operates  in  conjunction  with  the  SPS  executive  routines. 

IOLSQ  is  based  on  a  theory  employing  a  least- squares  fit  to  the 
computed  geocentric  position  vectors.  The  least  squares  procedure  does 
not  use  the  times  of  the  observations;  time  enters  into  the  formulation 
only  as  the  epoch  time,  i.e.,  the  time  of  the  first  observation.  As  the 
least  squares  method  is  not  iterative,  it  involves  only  linear  equations. 
The  result  is  an  initial  orbit  that  is  not  severely  affected  by  specific 
errors  due  to  a  few  pieces  of  data  This  is  in  contrast  to  an  exact  orbit 
based  on  only  two  position  vectors. 

Included  in  Section  Two  is  the  actual  formulation  used  in  IOLSQ 
to  provide  a  detailed  guide  to  the  program  operation.  Section  Three 
describes  the  input  and  output  requirements  while  the  flow  charts  are 
contained  in  Section  Four.  A  summary  of  the  experimentation  that  has  been 
conducted  with  the  program  to  evaluate  the  theory  is  included  in  Section 
Five. 

For  a  complete  discussion  of  the  differential  correction  portion 
of  IOLSQ,  see  Section  3  of  Reference  2  listed  at  the  end  of  this  document. 

A  discussion  of  the  ephemeris  part  of  differential  correction  can  be  found 
in  Reference  1. 
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SECTION  2 

FORMULATION 


Computation  for  the  Initial  Orbit  Least  Squares  Program  involves 
the  solution  of  a  number  of  linear  equations.  These  equations,  solved  in  a 
straight-forward  manner,  yield  the  classical  orbital  elements,  a,  e,  i,fi, 
and  CD.  Each  of  the  observations  must  initially  be  converted  to  position 
vectors,  in  the  inertial  system  thus  requiring  complete  observational 

fixes  of  range  and  two  angles  for  this  procedure.  The  formulation  follows: 

(1)  Determination  of  the  Geocentric  Position  Vector,  _r^,  and 
the  Unit  Position  Vector  U,~  : 

and  IL  are  calculated  for  each  observation  for 
i  =  1,2,3, ... ,N;  where  N  is  the  number  of 
acceptable  observations 

r.  =  p.L.  -  R.,  where  p.  is  the  observed  ranee. 

—  i  i—i  —i  i  ° 


r . 


The  relationship  between  the  satellite  position,  the 
dynamical  center,  and  the  observer  is  illustrated  in 
Figure  1. 
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FIGURE  1.  OBSERVATIONAL  GEOMETRY 


-3- 


The  position  of  the  Station,  R^ 


is  calculated  as  follows: 


x  =  fere]  cos  8 
Y  =  [dte]  stn  6 

z  =  z 

-\  /  2  2  2 
R  =  \JX  +  Y  +  Z 


where  Q,  the  local  sidereal  time,  is  computed  from 

0  =  0  +  X-, 

gr 


0gr  -  egr  +  D(0. 9856472)  +  F(360. 9856472) 

gr  =  Greenwich  sidereal  time  at  the  beginning  of  the  year, 
and  D  and  F  are  days  and  fractions  of  days,  respec¬ 
tively,  from  the  start  of  the  epoch  year  to  epoch. 

0  is  the  station  latitude  (positive  to  the  north) 


X  is  the  longitude  east 

hi 


- —  is  the  given  distance  of  the  station  from  the  polar  axis 

cos  y 

and  Z  is  the  polar  component  of  the  station's  geocentric  position. 
These  quantities  are  normally  acquired  in  referencing  the  SPS 
system  for  station  coordinates. 

The  unit  vector,  ,  is  computed  from  the  observed  quantities 
azimuth  (A)  and  elevation  (h) : 

L  =  -  cos  h  cos  A  sin  <h  cos  0  -  cos  h  sin  A  sin  0 
x 

+  sin  h  cos  0  cos  0 
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cos  h  cos  A  sin  0  sin  9  +  cos  h  sin  A  cos  9 


L  =  - 

y 


+  sin  h  cos  0  sin  9 


=  cos  h  cos  A  cos  0  +  sin  h  sin  0 

/  2  2  2 
such  that  L  =~\/  L  +  L  +  L 
Y  x  y  z 

(2)  Calculation  ofS  ,S  ,  S  ,S  ,  S  ,S  ,  the  Least  Squares 
_ xx'  xy’  yy’  xz’  yz’  zz’ _ _ 

Quantities  used  in  the  Solution  of  W  : 


S 

xx 


2 


N 


S 

xz 


S 

yz 


S 

yy 


N 


N 


I 


i  =  1 


U 


z  . 

l 


U 

yi 


u 

z . 

l 


s 

zz 
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(3)  Calculation  of  Wq,  the  Unit  Vector  Perpendicular  to  the  Orbit 
Plane : 

A 

W 

W  =  — t- 

W 


(See  Figure  2.) 


and 


A  ~\ 

/a  2 

A  2 

a  2 

W  =  \ 

J  W 

+  w 

+  W 

V  x 

y 

Z 

s 

s 

s  s 

A 

U  s 

_ 

xz  yy 

A 

W  = 

y 


A 

w 


s  s  -  s  ^ 

xx  yy  xy 


S  S  -  S  S 
xz  xy _ yz  xx 

S  S  -  S  2 

xx  yy  xy 


1.0 


(4)  Calculation  of  U  and  V  ,  where  V  is  the  Unit  Vector  in  the 
N  _ — o _ — Q _ O _ _ _ _ 

Orbit  Plane  Perpendicular  to  Uq: 


U  and  V  form  an  orthogonal  system  with  W 
— o  — o  ° 


Let 


Then  U 


n  = 

u.  •  w 

—1  — o 

A 

u 

U 

— 0 

__  — o 
”  A 

u 

o 

A 

u 

=  U 1  -  nW 

— 0 

-1  —o 

A 

Ia  2  A 

u 

=  U  +  U 

o 

/V  xo  y 

(subscript  1  denotes  the  first  observation) 


and 


V  =  w  x  U 
— o  — o  — o 
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FIGURE  2.  PROJECTION  OF  ORBIT-PLANE  ON  CELESTIAL  SPHERE  WITH  ORIENTATION 
UNIT  VECTORS  AND  ANGLES 
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Restricting  W  to  a  positive  value  will  generate  the  wrong 
zo 

solution  when  the  true  orbit  is  retrograde  and,  in  this  case, 
the  vector  V  will  be  calculated  in  the  wrong  direction.  To 
avoid  this  problem,  the  direction  of  motion  of  the  satellite 
is  determined  by  examining  a  pair  of  vectors,  and  U^+i> 

such  that  0  <  t.M  -  t.  <40  minutes.  If  W  is  selected 
l+l  l  —  z 

correctly,  then 

sign  of  (U^  x  =  sign  of 

where  j  denotes  the  vector  component  with  the  largest  magnitude. 

If  this  sign  equality  does  not  occur,  W  has  been  chosen 

o 

in  the  wrong  direction.  The  signs  of  the  components  of  Wq 
and  are  then  changed,  and  the  calculations  continue. 

(5)  Computation  of  cos(v^  -  vq)  and  sin  (v^  -  vq) : 


cos(v,  -  v  )  »  U.  •  U 
i  o  —i  —  o 


sin(v^ 


-  vo> 


U. 

—l 


V 

— o 


(6)  Formulation  of  S^,  S^,  S^j  S.^>  ^2’  S22’  S23’  ®3*  ^33’  t*ie  Least 
Squares  Quantities  used  in  the  Solution  of  p,  e  sin  vq  and  e  cos  vo: 


N. 

Sl  =  ^  sin  <Vi  '  V 
i=l 


SU  =  £  sln2  (vi  -  V 
i=l 
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I: 

[ 

I 

I 

I 

I 

1 

1 

I 

I 

I 

I 

I 

I 

■9 

I 

1 

L 

1 


12 


N 


sin  (v.  -  v  )  cos  (v.  -  v  ) 
1  o  i  o 


13 


N  . 

Y~  -  sin  (v . 

L  ri  1 

i=l 


V 

o 


) 


N 

S2  =  ^  cos  <Vi  -  vo) 

i=l 

N 

S22  =  Y_  C°s2  <Vi  '  Vo} 

l=i 


N 

S23  =  Y.  K  C°S  (Vi  ‘  V°} 

i=l 


N 


i=l 


N 


i=l 


2 
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e  cos  v  : 
o 


(7)  Calculation  of  p,  e  sin  vq> 


The  following  least  squares  equations  are  solved  to  obtain  p, 

e  sin  v  ,  e  cos  v  . 
o  o 


SL  =  (e  sin  vq)  -  (e  cos  vq)  S12  +  p  S13 

S«  =  (e  sin  v  )  S, ,  -  (e  cos  v  )  S00  +  p 
z  o  lz  o  ZZ  ZJ 


S3  =  (e  sin  v  )  S13  -  (e  cos  v  )  S23  +  p  S33 


(8)  Determination  of  e,  a,  e  cos  E  ,  e  sin  E  : 

o  o 

-  /  2  2 
e  =  y (e  cos  v  )  +  (e  sin  v  ) 

V  v  o  o 


e  cos  E 

o 


e  sin  E 

o 


e(e  +  cos  v  ) 


1  +  e  cos  v 

o 

Yl-e^  (e  sin  v  ^ 
o 


1  +  e  cos  v 


E 

o 


tan 


-1 


1-e^  sin  v 


e  +  cos  v 


o  J 


0  <  E  <2tt 
—  o 
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(9)  Calculation  of  h,  the  Angular  Momentum  Vector  Directed  Along  W  : 


h  =  nTp  wo 


(!°)  Calculation  of  12  and  i,  the  Longitude  of  the  Ascending  Node  and 
the  Inclination  of  the  Orbit  Plane: 

Since  the  three  equatorial  component  equations  of  W  are 

— o 


W  =  sin  (1  sin  1 
x 

o 


cos  12  sin  i 


W  =  cos  i 
z 

o 

then  12  =  tan 


and  i  =  tan 


0  <  12  <  2tt 

0  <  i  <  f 


(11)  Calculation  of  a,  the  Vector  Directed  to  the  Perifocus,  and  Its 

Components,  a  „  and  a  •  ~  '  ’  “ 

_  xN _ 

a  ^  is  the  component  of  a  along  N,  the  unit  vector  directed  to 
ascending  node. 

is  the  component  of  a  along  M,  the  unit  vector  in  the  orbit 
plane  perpendicular  to  N,  in  the  direction  of  motion. 
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(12)  Calculation  of  10,  the  Argument  of  Perigee: 


cjd  =  tan 


0<  CD  <  2  7T 


(13)  Calculation  of  fi,  co  ,  Pa>  and  : 

Q  is  the  rate  of  change  of  the  longitude  of  ascending  node. 

00  is  the  rate  of  change  of  the  argument  of  perigee. 

Pa  is  the  anomalistic  period. 


H 

q 


is  the  altitude  of  perigee  in  statute  miles. 


k  W 

•  e  z 

°  =  1/2  "2 -  (1.62329  x  10'J) 

a  p 


O) 


- Tjl - 2  <L- 62329  x  10’3)  (5  Wz  2  -  1) 

2  a  p  o 
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(14) 


_  2TT_  3/2 

a  k  a 
e 

Hq  =  [a  (1-e)  -  1]  (3963.199) 

Calculation  of  the  Mean  Longitude  at  epoch,  L  : 

L  =  +  0  +  CD  +  E  -  e  sin  E  if  i<  90° 

o  —  o  o  — 

if  i  >  90° 


use  +  Q 
use  -  Q 
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SECTION  3 


PROGRAM  OPERATIONS 


This  program  is  executed  in  conjunction  with  the  B-2  Semi- 
Automatic  Program  System  (SPS)  job  schedule  mode  of  operation.  The  normal 
operating  procedures  required  for  this  mode  of  program  execution  are 
described  in  Reference  3. 

Observations  are  read  into  memory  from  cards.  Sensor  coordinates 
are  stored  in  core  from  cards  or  from  the  system  SEAI  file.  The  position 
vectors  of  the  satellite  are  then  calculated  and  stored  in  memory. 


3 . 1  INPUT 


IOLSQ  has  two  input  options: 

0  :  Observation  Cards. 

1  :  Observation  and  Sensor  Cards. 

Required  input  data  consists  of  observation  types  2  cr  3  and/or 
the  Sensor  Cards.*  These  data  cards  must  be  used  in  conjunction  with  the 
usual  SPS  schedule  tape  mode  control  cards. 


Observation  Types: 


2  -  contains  elevation,  azimuth  and  range 

3  -  contains  elevation,  azimuth,  range  and  range-rate. 
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The  tape  setup  required  for  input  is  arranged  as  follows: 


TAPE 

1 

2 

4 

11 

The  Schedule  Tape 
punched  cards: 

ID  Card 
JOB  Card 
REM  Card 
S PS JOB  Card 
Sensor  Cards 
Observation  Cards 
END  Card 
END  Card 
Blank  Card . 


DESCRIPTION 
SPS  Master  Tape 
Schedule  Tape 
SEAI  Tape 
Output  Tape 

used  for  input  is  generated 
(70  SC  TP) 

( IOLSQ) 

-  optional 

(END  OF  JOB) 

(ENDSCHED) 


from  the  following 


The  items  shown  in  parentheses  above  are  the  code  names  to  be  used  in  the 
respective  cards.  These  cards,  along  with  the  JOB,  REM,  Sensor  and  Obser  - 
vation  Cards  are  described  in  Section  5.4  of  Reference  2. 


3 . 2  OUTPUT 

All  printed  or  punched  output  may  be  obtained  from  logical  tape 
unit  11,  using  data  select  1  for  printed  output  and  data  select  2  for 
punched  card  output  . 
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The  printed  output  of  IOLSQ  begins  with  a  heading  and  the  date. 
This  is  followed  by  the  N,  M  Elements  and  the  Classical  Elements  at  the 
time  of  the  first  observation  (epoch  time).  The  differential  correction 
output  follows,  showing  the  observation  residuals  and  correction  data  at 
each  iteration  followed  by  the  final  correction  summary.  Printout  for 
the  differential  correction  portion  is  exactly  like  that  of  all  SPS 
differential  correction  packages.  (See  Reference  2.).  In  this  standard 
differential  correction  output,  the  corrected  elements  are  displayed, 
plus  the  corrected  elements  updated  to  the  final  revolution.  For  this 
program,  the  results  of  the  two  sets  of  new  elements  will  be  the  same. 

See  Figure  3,  page  19 . 

The  punched  card  output  consists  of  two  standard  SPS  7-Card 
Element  sets.  System  routines  were  used  to  obtain  the  punched  output 
so  that  the  two  sets  are  duplicates.  This  is  because  the  corrected  set 
and  the  updated  set  are  the  same  in  IOLSQ. 


-16- 


3.3 


PROGRAM  NOTES 


Comments  may  be  printed  indicating  the  operation  of  the  IOLSQ 
Program.  Any  of  the  comments  from  the  following  list  may  be  produced 
during  the  operation  of  the  initial  orbit  portion  of  IOLSQ: 

T  FFT 

(1)  "EXPONENT  OVERFLOW  AT  THE  HALF  OF  LOC.  XXXXX" 

RIGHT 

T  FFT 

(2)  "SUBROUTINE  ERROR  AT  THE  “1“  HALF  OF  LOC.  XXXXX" 

RIGHT 

(3)  "NO  COORDS  SPECIFIED  FOR  SENSOR  NO.  XXXX" 

This  comment  indicates  that  the  OBSGET  subroutine 
could  not  find  from  either  the  SEAI  Tape  or  from 
the  Sensor  Cards  the  coordinates  corresponding  to 
the  sensor  number  on  one  of  the  Observation  Cards. 

(4)  "THERE  ARE  XXX  OBSERVATIONS  WHICH  HAVE  NO  RANGE" 

Only  observations  types  2  and  3  can  be  used  by  the 
initial  orbit  part  of  IOLSQ.  This  comment  indicates 
the  number  of  input  observations  that  were  rejected 
during  the  initial  orbit  determination. 

(5)  "NUMBER  OF  GOOD  OBS .  IS  LESS  THAN  5" 

The  least-squares  theory  in  IOLSQ  requires  that 
there  be  a  minimum  of  5  acceptable  observations 
in  order  that  the  calculations  can  proceed. 

(6)  "OBSERVATION  BUFFER  FULL,  ONLY  THE  FIRST  150  TYPE 
2  or  3  INPUT  OBS.  WILL  BE  USED" 

This  comment  indicates  that  the  maximum  of  150 
acceptable  observations  was  exceeded.  Therefore, 
calculations  will  proceed  using  only  the  first 
150  type  2  and  3  observations  . 

(7)  "THERE  IS  NO  PAIR  OF  OBSERVATIONS  WHICH  ARE  SEPAR¬ 
ATED  IN  TIME  BY  MORE  THAN  1  MIN.  AND  LESS  THAN  40 
MIN." 
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To  determine  the  direction  of  motion  of  the 
satellite,  the  program  attempts  to  examine  two 
observations  which  are  greater  than  1  minute 
but  less  than  40  minutes  apart.  The  above 
comment  indicates  that  no  such  pair  was  found 
and  the  job  is  terminated. 

The  error  comments  that  may  result  from  the  differential  correc¬ 
tion  portion  of  the  program  are  fully  described  in  Section  3,  pages  3-84 
and  3-85,  of  Reference  2. 


! 
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FIGURE  3.  IOLSQ  OUTPUT  FORMAT 


SECTION  4 


FLOW  DIAGRAMS 


The  following  illustrations  indicate  the  computational  procedures 
of  IOLSQ  in  diagrammatic  form.  Diamonds  indicate  decisions;  arrows  indi¬ 
cate  the  direction  of  logical  flow;  ovals  indicate  subroutines,  and  rectangles 
indicate  mathematical  or  bookkeeping  operations.  The  small  circles  indicate 
starts  and  exits  and  also  serve  as  connectors  from  one  page  to  the  next.  The 
oval  at  the  entrance  to  a  subroutine  contains  the  name  of  that  particular 
subroutine.  A  rectangle  connected  to  a  subroutine  oval  contains  a  brief 
description  of  that  particular  subroutine.  Dotted  lines  from  a  subroutine 
oval  indicate  possible  exits. 
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FIGURE  4.  IOLSQ  GENERAL  FLOW  DIAGRAM 
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FIGURE  5.  IOLSQ  DETAILED  FLOW  DIAGRAM  (1  OF  3) 
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EACH  OBS.  AND 
STORE  IN  TBLOC 
BUFFER 


COMPUTE  a,  a  ,  a  CD: 

—  xN  yN 

a  =  U  e  cos  v  -  V  e  sin  v 
—  — o  o  — o  o 


a  W  -  a  W 
y  x  x  y 


FORM  THE  FOLLOWING  3  LEAST- SQUARES 
EQUATIONS  FROM  TBLOCK: 

S1  =  P  Sn+  [e  sin  vq  ]  Su-  [e  cos  vq  ]  S12 

s 2  ^  P  s23+  Ce  sln  vo  ]  S12‘  COS  Vo  1  S 22 

s 3  -  P  s33+  Si°  Vo  ]  S 13~  [e  COS  vo  ]S23 


-W  [a  W  +  a  H  ]  ; - tt 

2  x  *  r  -4  a  Vi-w2 

/  2  z  z 

v  i  -  w  * 

z 


00  =  tan 


■■ft] 


FIGURE  5.  IOLSQ  DETAILED  FLOW  DIAGRAM  (2  OF  3> 
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.  IOLSQ  DETAILED  FLOW  DIAGRAM  (3  OF  3) 
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FIGURE  5 


SECTION  5 

PROGRAM  EXPERIMENTATION 


5.1  SUMMARY  OF  TESTING 

In  this  section  a  summary  of  twenty  cases  that  were  run  with 
the  IOLSQ  Program  is  presented  to  illustrate  the  capabilities  of  the 
initial  orbit  determination. 

Observations  were  chosen  from  four  different  satellites,  each 
satellite  representing  a  different  type  of  orbit.  The  following  table 
summarizes  the  different  orbits  used. 

TABLE  I 

SUMMARY  OF  TEST  ORBITS 
Satellite 


imber 

Drag 

Eccentric ity 

Oil 

Low 

Moderate 

070 

Moderate 

Low 

116 

Low 

Low 

213 

High 

Low 

Five  runs  were  made  for  each  satellite  with  each  successive 
run  having  more  observations  than  the  previous  run  and  covering  a  longer 
time  span.  The  time  span  and  the  number  of  observations  used  are 
included  in  the  following  tables  along  with  the  initial  and  corrected 
elements.  The  test  cases  are  presented  here  to  illustrate  the  accuracy 
of  the  initial  orbit  portion  of  IOLSQ  rather  than  the  differential  correc¬ 
tion  portion.  Each  initial  orbit  is  compared  with  a  reference  orbit  which 


-25- 


was  obtained  from  the  SGPDC  Program.  (See  Reference  2.)  The  reference 
orbit  was  determined  using  file  elements  and  the  same  observations  as 
the  initial  orbit  case  with  the  largest  number  of  observations. 


5 . 2  RESULTS 


Experimentation  indicates  that  IOLSQ  is  quite  adequate  for 
determining  an  initial  orbit  when  a  group  of  associated  observations 
is  given.  The  initial  orbital  parameters  are  entered  as  input  to  the 
differential  correction  scheme,  and  the  program  appears  to  yield  suffi¬ 
ciently  accurate  elements  for  this  purpose. 

In  the  early  phase  of  the  experimentation,  the  initial 
elements  were  not  always  sufficiently  accurate  when  corrected  over  a 
long  time  span  of  observations.  This  led  to  divergence  in  the  correc¬ 
tion  process.  To  mitigate  this  problem,  a  special  period  correction 
was  included  in  IOLSQ  to  improve  the  semi-major  axis.  This  correction 
is  executed  by  fitting  a  straight  line  in  the  least-squares  sense  to 
the  time  residuals.  The  correction  to  the  semi-major  axis  is  then 
computed  from  the  slope  of  this  line. 

When  the  time  span  covered  by  the  observations  is  long,  the 
time  residuals  may  exhibit  a  "saw-tooth"  effect.  The  time  residual 
is  similar  to  the  angle  in  the  orbit  plane  between  the  computed 
position  and  the  observed  position.  The  time  residual  is  thus  computed 
from  the  angle  between  two  vectors  and  necessarily  lies  between  -  7 T  and 
+  7T  .  Therefore,  when  the  time  residuals  exceed  one-half  the  period, 
they  will  change  sign  and  exhibit  the  saw-tooth  effect.  The  sketch 
below  illustrates  A  t  as  a  function  of  time: 


* 

i 
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To  fit,  in  a  least-squares  sense,  a  line  to  the  time  residuals, 
it  becomes  necessary  to  eliminate  the  "saw-tooth",  i.e.s  A  t  must  be  made 
a  linear  function  of  time.  This  is  accomplished  by  adding  or  subtracting 
the  value  of  the  period  when  two  successive  time  residuals  differ  by  more 
than  one-half  the  period.  This  special  period  correction  is  invalid  in 
the  case  where  two  successive  time  residuals  are  far  enough  apart  that 
there  exists  more  than  one  change  of  sign  between  them 

IOLSQ  differs  from  other  initial  orbit  programs  in  that  the 
semi-major  axis  is  determined  from  geometrical  rather  than  dynamical 
considerations.  Thus,  it  is  to  be  expected  that  there  will  be  an  error 
in  the  semi-major  axis  when  the  points  in  space,  as  computed  from  the 
observations,  do  not  uniquely  define  an  ellipse  The  experimentation 
has  shown  that  the  semi-major  axis  is,  indeed,  inadequately  determined 
when  the  observations  are  concentrated  at  one  point  in  the  orbit  but 
becomes  increasingly  better  as  more  observations  are  added.  This  effect 
is  illustrated  in  the  3-hour  correction  of  Satellite  Oil  (see  Table  II). 

The  initial  orbit  formulation  assumes  that  the  longitude  of 
the  ascending  node,  (2  ,  is  constant.  It  can  therefore,  be  expected  that 
due  to  the  precession  of  the  node,  £)  will  be  poorly  determined  when  the 
span  of  the  observations  is  large.  This  limitation  of  the  program  is 
illustrated  by  Satellites  070  and  213  (see  Tables  III  and  V)  where  fi 
becomes  increasingly  worse  and  is  in  error  bv  as  much  as  16°  when  the 
observation  span  is  16  days  and  15  hours . 

Another  limitation  of  the  program  is  illustrated  by  Satellites 
070,  116  (see  Table  IV),  and  213  The  argument  of  perigee, CD ,  is  poorly 
defined  when  the  eccentricity  is  low  The  large  errors  in  co  are  due  to 
the  low  eccentricity  and  should  be  corrected  by  the  differential  correc¬ 
tion  For  low  eccentricities,  however,  a  large  error  in  CD  does  not 
necessarily  indicate  an  inaccurate  orbit. 

The  eccentricity  and  the  argument  of  perigee  are  poorly 
defined  when  the  observations  represent  a  single  point  in  space  In 
this  case,  the  program  cannot  uniquely  determine  an  ellipse  to  fit  the 
observations.  This  situation  is  illustrated  by  Satellite  Oil,  Case  5 
The  observations  are  concentrated  at  a  single  point  in  space,  and  the 
differential  correction  could  not  adequately  determine  an  orbit  except 
when  the  correction  of  the  eccentricity  and  argument  of  perigee  was 
inhibited 


The  effect  of  high  drag  is  illustrated  by  Satellite  213 
The  computation  of  drag  is  not  included  in  IOLSQ.,  and  consequently 
the  differential  correction  cannot  always  converge  to  an  orbit  that 
closely  fits  the  observations .  Cases  3,  4  and  5  could  converge  to  an 
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orbit  only  when  the  correction  of  the  eccentricity  was  inhibited.  When 
the  eccentricity  was  corrected,  the  differential  correction  began  to 
diverge  steadily.  These  cases  represent  the  best  results  that  can  be 
obtained  using  the  usual  "no-drag"  assumption. 

The  effect  of  bad  observations  is  illustrated  by  the  cases 
run  on  Satellite  070.  Cases  3,  4  and  5  contain  a  group  of  observations 
that  represent  a  single  point  in  space.  This  point  is  net  compatible 
with  the  rest  of  the  orbit,  and  the  differential  correction  could  not 
determine  an  orbit  to  fit  all  of  the  observations  except  when  the 
eccentricity  was  not  corrected.  Cases  4  and  5  show  that  the  effect 
of  these  bad  observations  is  minimized  when  observations  are  added 
that  more  clearly  define  the  orbit. 

As  a  result  of  the  testing  (and  in  view  of  the  foregoing 
comments) ,  it  is  concluded  that  IOLSQ  is  best  suited  for  orbit  deter¬ 
minations  using  groups  of  observations  that  cover  a  time  span  on  the 
order  of  four  or  five  days.  Enough  observations  must  be  available  for 
input  to  enable  the  program  to  geometrically  determine  the  semi-major 
axis.  The  time  span  must  be  small  enough  so  that  the  precession  of  fi 
is  not  a  dominating  factor.  The  advantage  of  IOLSQ  is  that  its  accuracy 
is  not  severely  affected  by  errors  in  a  few  observations  as  would  occur 
for  other  initial  orbit  programs*  based  on  one-  or  two-position  fixes. 


5.3  EXAMPLE  OF  PROGRAM  OPERATION 

Portions  of  an  IOLSQ  test  case  are  presented  on  the  following 
pages  to  illustrate  the  combined  operation  of  the  initial  orbit  deter¬ 
mination  and  the  differential  correction  procedure.  The  test  case  is 
reproduced  in  the  actual  computer  output  format  for  Satellite  116, 

Case  3.  (See  Figure  7.)  Case  5  printed  at  the  end  of  the  output 
indicates  the  fifth  differential  correction  iteration  for  Case  3.  Due 
to  space  limitations,  the  test  case  output  is  abbreviated  to  show  the 
initial  elements,  the  last  pass  residuals,  and  the  final  element  summary. 
All  input  data  for  this  case  are  shown  in  the  reproduced  80-80  listing 
of  the  input  deck  (see  Figure  6). 


*Such  as  Initial  Orbit  Radar  Fix  ( IORF)  Program,  which  computes  the 
initial  orbit  from  2  or  3  radar  fixes,  and  Initial  Orbit  Radar  Track 
(I0HG)  Program,  which  calculates  the  initial  orbit  from  a  radar  track. 
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★All  observations  except  one  within  3  min.  of  one  another,  hence  the  poor  semi-major  axis. 
*e  and  CDnot  corrected  in  the  differential  correction  process. 
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FIGURE  6.  SAMPLE  INPUT  DATA 
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FIGURE  7.  SAMPLE  OUTPUT  DATA  -  SATELLITE  116,  CASE  3  ( 1  OF  5) 
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